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AUTOMATED  SYSTEM  FOR  ANALYSIS  OF 

RUNOFF  HYDROGRAPHS 

\     I         $~- 
~7a  ^  WjD.  L. iBrakensiek \, 

INTRODUCTION 

Hydrologists  often  utilize  dimensionless  hydrographs  to  predict  simple  hydrographs  (single 
peak  with  the  recession  not  sustained  by  rainfall  or  return  flow).  The  dimensionless  hydrograph 
as  introduced  in  this  paper  has  a  discharge  scale  of  (q/qp)  and  a  time  scale  of  (t/tp)  where 

q  =  runoff  rate 

qp  =  peak  runoff  rate 

t  =  time 

tp  =  time  to  peak  rate 

By  using  the  assumptions  of  superposition  and  linearity,  hydrologists  develop  complex  hydrographs 
from  a  summation  of  simple  hydrographs. 

Hydrographs  for  simple  runoff  events  have  been  noted  to  be  similar  in  shape  to  skewed  sta- 
tistical frequency  functions.  Gray  2  and  Reich  utilized  this  justification  for  fitting  a  gamma 
function  to  a  number  of  observed  hydrographs.  Relating  the  fitted  function  parameters  to  hydro- 
logic  and  geomorphic  variables  makes  it  possible  to  predict  hydrographs  for  ungaged  watersheds. 

This  study  considers  a  function  similar  to  the  gamma  function,  except  that  it  has  only  one 
parameter.  Even  though  the  function  is  a  one-parameter  function,  it  is  used  in  a  two-parameter 
fashion,  that  is,  two  values  of  the  parameter  are  estimated,  one  for  the  rising  and  one  for  the 
falling  side  of  the  hydrograph.  A  computer  program  was  written  to  automate  the  fitting  procedure. 
A  program  listing  and  example  are  presented  in  the  appendix. 

The  function  was  fitted  to  100  hydrographs  and  the  quality  of  fit  was  evaluated  from  the  de- 
viation of  the  estimated  hydrographs  from  the  observed  hydrographs. 


HYDROGRAPH  FUNCTION 


As  in  any  instance  of  curve  fitting,  only  empirical  justification  can  be  presented  for  the  fitting 
function.  The  function  developed  for  this  study  is  written  as 

q/qp  =  EXP    ^-n  (In  (t/tp) +  2  (l/v/t/tp"  -1)))  [1] 
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or,  if  we  define  a  transformed  time  scale  as 
T=  In  (t/tp)+  2  (lVt/t"  -1) 

then  equation  1  is  written  as 

q/qp  =  EXP   (-nT)  [2] 

where  n  is  the  parameter  to  be  estimated.  Equation  2  is  now  similar  to  the  usual  exponential 
recession  function.  The  function  1  is  fitted  to  the  rising  and  falling  side  of  the  dimensionless 
hydrograph.  Hence,  a  two-parameter  system  is  fitted.  Of  course,  estimates  of  qp  and  tp  are 
required  to  recapture  the  natural  hydrograph. 


FUNCTION  FITTING 

The  Fortran  II  source  program  is  listed  in  Appendix  I.  It  was  written  for  an  IBM  1620  com- 
puter (40  K  memory  with  indirect  addressing).  The  program  performs  essentially  four  types  of 
calculations. 

(1)  The  tabulated  time  scale  (civil  or  military  time)  of  the  natural  hydrograph  is  converted  to 
decimal  hours.  The  decimal  hour  scale  is  then  converted  to  an  accumulated  hourly  time 
scale. 

(2)  A  search  is  made  for  the  peak  rate  and  time-to-peak  rate  of  the  natural  hydrograph.  The 
time  and  discharge  scales  of  the  natural  hydrograph  are  then  converted  to  dimensionless 
scales. 

(3)  The  dimensionless  hydrograph  scales  are  transformed  in  accordance  with  the  fitting 
function: 

(a)  Discharge  scale  =  In  (q/qp) 

(b)  Time  scale  =  In  (t/tp)  +  2  U/Vt/tp  -1) 

(4)  The  rising  and  falling  side  parameters  are  estimated  by  the  usual  least  squares  formula 
applied  to 

ln(q/qp)=  -nT. 

It  was  apparent  from  manual  fitting  trials  that  a  closer  fit  is  more  important  in  high-flow 
parts  of  the  hydrograph  than  in  the  extreme  tail  parts  of  the  hydrograph.  To  incorporate  this  into 
the  fitting  program,  a  cutoff  constant  was  utilized.  By  varying  its  value,  the  portion  of  the  hydro- 
graph  contributing  to  estimation  of  the  parameter  is  controlled.  Values  of  this  constant  were 
from  3.0  in  decrements  of  0.5  down  to  0.5.  In  terms  of  cutoff  of  the  (q/qp)  -  scale,  the  equivalence 
is,  for  example, 


C  =  1. 

,   (q/qp)  =  0.367 

C  =  2. 

,   (q/qp)  =     .135 

C  =  3. 

,   (q/qp)  =    .049 

Thus,  using  a  cutoff  of  2.0,  the  function  is  fitted  to  values  of  (q/qp)  between  0.135  and  1.00.  The 
program  simply  selects  the  closest  fit  on  the  basis  of  minimum  sums  of  squares  of  observed 
minus  calculated  from  the  sequence  of  C-values. 


4 
Mention  of  a  trade  name  does  not  constitute  a  guaranty  or  warranty  of  the  product  named  and  does  not  signify 

that  this  product  is  approved  to  the  exclusion  of  other  comparable  products. 


The  final  output  of  the  program  provides  the  best  fitting  hydrograph,  the  observed  hydro- 
graph,  time  to  peak,  peak  rate,  and  the  rising  side  and  falling  side  parameter  values. 

In  Appendix  II  the  data  input  formats  are  given.  For  each  hydrograph,  one  parameter  card 
is  required  (C  is  usually  valued  initially  as  3.).  One  data  card  is  required  for  each  line  of  data. 

Appendix  III  presents  a  sample  set  of  input  data  and  the  corresponding  program  output.  Ex- 
perience to  date  indicates  that  20  to  30  hydrographs  can  be  fitted  per  hour,  depending  on  the 
number  of  entries  per  hydrograph. 


APPLICATION 


An  initial  study  was  undertaken  to  assess  the  fitting  value  of  this  hydrograph  function.  The 
study  utilized  100  published  hydrographs  found  in  "Selected  Runoff  Events  from  Small  Agricul- 
tural Watersheds  in  the  United  States"  5  and  "Hydrologic  Data  for  Experimental  Agricultural 
Watersheds  in  the  UnitedStates,  1956-59".6  Table  1  presents  the  watershed  and  storm  identifica- 
tion together  with  the  estimated  values  of  the  parameters. 


Table  1. --Listing  of  watersheds  and  storms   studied  and 

fitted  parameters. 


Location  and  watershed 

Storm 
Month- day-year 

n  Rise 

n  Fall 

Bentonville,  Ark. : 

W-5 

5-25-39 

-5.4313 

-7.9732 

W-5 

6-9-52 

-34.3306 

-20.7748 

Safford,  Ariz. : 

W-I 

7-19-57 

-1.4565 

-10.2684 

W-II 

8-20-56 

-3.0521 

-7.8796 

W-II 

7-16-59 

-4.3990 

-20.6930 

w-v 

8-28-57 

- . 6709 

-3.3382 

W-V 

8-30-57 

-.4415 

-2.6094 

Tombstone,  Ariz. : 

W-2 

7-20-59 

-6.2263 

-1.6997 

Colorado  Springs,  Colo. : 

W-IV 

8-13-45 

-33.5999 

-22.2823 

Americus,  Ga. : 

W-IV. 

8-19-42 
1-17-43 

-39.5155 
-11.8268 

-14.8924 

W-IV 

-7.0215 

U.S.  Agricultural  Research  Service.  Selected  Runoff  Events  for  Small  Agricultural  Watersheds  in  the  United 
States.   U.S.  Dept.  Agr.,  Agr.  Res.  Serv.,  Soil  and  Water  Conservation  Research  Division,  374  pp.   January  1960. 

Hobbs,  Harold  W.,  comp.  Hydrologic  Data  for  Experimental  Agricultural  Watersheds  in  the  United  States, 
1956-59.    U.S.   Dept.   Agr.    Misc.  Pub.  No.  945,    672  pp.    November  1963. 


Table  1. --Listing  of  watersheds  and  storms  studied  and 
fitted  parameters --Continued. 


Location  and  watershed 


Storm 
Month-day-year 


n  Rise 


n  Fall 


Edwards ville,  111. : 

W-l ■ 

W-4 

Iowa   City,    Iowa: 
Ralston   Creek 

College  Park,   Md. : 

W-l 

W-l 

W-2 

W-6 

W-6 

W-7 

W-7 

Oxford,    Miss.  : 

W-4 

W-4 

W-5 

W-5 

W-5 

W-10 

W-28 

W-28 

W-35 

WC-1 

WC-2 

WC-3 

Hastings,    Nebr. : 

W-3 

W-3 

W-5 

1-H 

1-H 

2-H 

4-H 

Albuque  r que ,    N .    Mex . 

W-I 

W-I 

W-I 

W-I 

W-II 

W-III 


5-27-38 
5-27-38 


7-18-56 


7-22-45 
6-15-54 
7-22-45 
8-27-43 
7-22-45 
8-10-42 
7-22-45 


4-3-58 

9-9-59 

4-3-58 

6-10-59 

6-11-59 

5-22-59 

7-22-58 

9-9-59 

5-22-59 

5-26-59 

5-26-59 

5-26-59 


5-1-57 
6-15-57 

7-3-59 
6-16-57 
6-12-58 
6-12-58 

5-4-59 


9-8-47 
8-19-56 

8-9-57 
8-14-58 
8-24-57 
8-19-56 


-83.8372 
-90.0152 


-90.6247 


-46.3302 
■11.1429 
-9.3360 
-5.2201 
-15.0175 
■31.7859 
-13.4699 


-21.2269 

-3.1856 

-8.0191 

-10.1031 

-4.3152 

■113.8592 

-25.3347 

-4.1932 

-19.7394 

-16.9993 

-10.0764 

-26.7372 


-45.1365 

■13.6296 
-9.8426 

•31.5879 
-5.9608 
-8.8988 

■11.9878 


-20.4245 
■18.2309 

-3.4505 
■32.3868 
■25.8804 

-5.6182 


-40.3424 
-21.7968 


-10.8163 


■32.6977 
■18.2174 
■14.5556 
■16.5945 
■14.3344 
■27.0364 
■31.3912 


-8.2025 

-10.1465 

-7.8754 

-9.7004 

-6.4076 

■27.1027 

-15.0764 

-5.9237 

■17.0422 

-30.4440 

-19.8004 

-32.9694 


■12.8659 
■18.5796 
■12.9281 
-9.4677 
•10.9489 
•14.8207 
-5.7904 


•13.0889 

-4.6776 

-7.7425 

-46.9678 

-11.8369 

-8.4936 


Table  1. --Listing  of  watersheds  and  storms  studied  and 
fitted  parameters --Continued. 


Location  and  watershed 


Storm 
Month- day-year 


n  Rise 


n  Fall 


Mexican  Springs,  N.  Mex. 

W-2 

W-ll 


High  Point,  N. C. : 

West  Fork  of  Deep  River. 


Coshocton,  Ohio : 

131 

127 

110 

Ill 

187 

169 

183 , 

196 , 

94 , 

Hamilton,  Ohio: 

W-l , 

W-l , 

Cherokee,  Okla. : 

W-l , 

W-i , 

W-l , 

W-2 

W-2 , 

W-3 

W-4 , 

W-4 

Guthrie,  Okla. : 

W-VI 

W-VI 

Stillwater,  Okla. 

W-l 

W-l 

W-3 

Riesel,    Tex. : 

W-l 

W-l 


9-5-40 
8-26-39 


5-2-39 

12-26-38 

2-26-39 


6-12-57 

6-12-57 
6-12-57 
6-12-57 
6-12-57 
6-12-57 
8-16-47 
8-16-47 
6-12-47 


5-17-43 
7-7-43 


6-24-58 
9-25-59 
5-28-60 
6-24-58 
10-13-59 
5-28-60 
6-24-58 
5-28-60 


9-8-42 
6-26-45 


4-18-57 
6-27-57 
6-27-57 


6-10-41 
3-26-46 


-23.3866 
-1.8471 


-14.0067 
-6.4381 
-3.2270 


-2.7887 

-184.6239 

-3.0373 

-50.3599 

-59.3227 

-83.1074 

-15.9781 

-313.7573 

-44.8717 


■10.0951 
-37.8762 


-29.3441 

-7.1256 

•14.0688 

■36.5202 

■29.5597 

-8.3606 

-5.3942 

-4.6167 


-24.4832 
-37.9473 


-93.7714 

-38.4060 

-4.1090 


-33.4406 
■195.0137 


-39.2123 
-6.9189 


-4.4998 

-20.8938 

-5.6469 


-10.3692 
-19.9397 
-25.3157 
-39.2510 
-11.3100 
-43.8798 
-15.1650 
-97.4145 
-32.7548 


-3.3781 
■10.8566 


-8.7506 
-13.5323 
-42.1781 
-21.4625 
-20.0797 
-23.3567 

-6.3858 
-21.8304 


-7.7950 
-8.0815 


■341.2761 
•153.5767 
-12.6162 


-19.3325 

-64.3905 


Table  1. --Listing  of  watersheds  and  storms  studied  and 
fitted  parameters-- Continued. 


Location  and  watershed 


Storm 
Month- day-year 


n  Rise 


n  Fall 


Riesel,  Tex.— Continued: 

W-l ' , 

W-l 

W-2 , 

W-6 , 

W-6 , 

W-10 , 

W-10 , 

Y-6 

Y-6 

Y-6 

SW-12 , 

SW-17 

SW-17 

Staunton,   Va. : 
Bell  Creek 
W-l 

Fennimore,   Wis. : 

W-l 

W-2 

W-3 

W-3 

W-3 

W-3 

LaCrosse,   Wis. : 

CW 

CW 

CW 


4-27-49 
4-24-57 
6-23-59 
4-24-57 
6-23-59 
4-24-57 
6-23-59 
4-24-57 
6-4-57 
6-23-59 
6-23-59 
3-31-57 
4-24-57 


9-8-48 


8-5-51 
8-5-51 
8-12-43 
7-11-44 
6-28-45 
6-24-49 


8-16-40 
6-29-41 
9-15-41 


-29.2257 
■262.7649 
•310.2265 

-7.9467 
-63.8263 
-10.4061 
-83.3682 
-17.9326 

-6.0614 
-33.0627 
-38.7213 
-24.2796 

-7.9012 


-22.5718 


-84.5015 

-12.0619 

-9.7027 

-455.6838 

-45.4637 

-1.5383 


-14.1442 

-3.0925 

-21.7221 


-33.9124 

■117.3924 

-81.1602 

-13.9224 

-22.4897 

-15.1351 

-53.8528 

-6.6444 

-5.8565 

-27.4650 

-53.1431 

-6.5723 

-17.3511 


■10.8254 


■99.8869 

■25.4745 

-6.2509 

-10.9080 

-9.3694 

-1.5163 


-10.9085 

-5.1205 

■27.1035 


The  fitting  ability  of  the  hydrograph  function  was  evaluated  from  each  of  the  storms  listed  in 
table  1.  As  mentioned,  for  each  estimated  parameter  there  is  a  sum  of  squares  of  deviation  for 
the  difference  between  the  observed  and  estimated  hydrograph  points.  This  was  printed  out  for 
this  evaluation.  All  of  these  were  grouped  for  the  rising  and  falling  parameters  according  to 
magnitude.  In  table  2  the  total  frequency  distribution  of  the  quality  of  fit  is  presented.  Reference 
is  made  to  figures  1  and  2  for  quantitative  measures  of  the  quality  of  fit. 

EXCELLENT  Figure  1,  rising  and  falling  side,  is  an  excellent  fit. 

GOOD  Figure  2,  rising  side,  is  at  the  lower  end  of  the  good  fit  scale. 

FAIR  TO  POOR  Figure  2,  falling  side,  is  a  poor  fit. 

Such  hydrographs  as  figure  2  should  really  not  be  fitted,  since  the  rising  side  is  character- 
istic of  an  initial  runoff  period  followed  by  the  main  runoff  period,  and  the  recession  side  is 
somewhat  sustained  by  a  late  occurrence  of  rainfall  or  return  flow.  Some  evidence  of  better 
fitting  ability  on  the  falling  limb  of  the  hydrograph  is  shown  in  table  2. 
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Figure  1. — Hydrograph  fit. 
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Figure  2. — Hydrograph  fit. 
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Table  2. --Evaluation  of  the  fitting  system. 


Percent  of  cases 

Quality  of  fit 

Rising  side 

Falling  side 

28.0                  30.0 

Good2. 

46.0                  53.0 

26.0                 17. n 

-'■Both  sides  of  Figure  1. 

2Rising  side  of  Figure  2  is  at  the  lower  end  of  the  good  fit  scale, 

3FaIling  side  of  Figure  2  is  a  poor  fit. 


CONCLUSIONS 


To  date  no  particular  difficulties  have  been  encountered  in  the  application  of  the  computer 
program. 

The  function  presented  in  this  study  for  fitting  dimensionless  hydrographs  has  performed 
satisfactorily  when  applied  to  strictly  simple  hydrographs.  The  poor  fits  often  resulted  from  using 
complex  hydrographs. 

The  success  of  this  study  has  made  it  feasible  to  examine  a  much  larger  population  of  simple 
hydrographs.  Development  of  parameter  prediction  equations  can  then  proceed.  Once  these  are 
developed,  it  would  be  possible  to  estimate  simple  hydrographs  much  in  the  same  fashion  as  unit 
hydrographs  are  estimated. 


SYMBOLS  USED  IN  THE  PROGRAM 


Symbol 
TIME  (I) 
TIM 

TIME1  (I) 
TPEAK 
PEAK 
R 
T 

DISCH1  (I) 
TIME2  (I) 
TSCALE  (I) 


Description 


Time  in  decimal  hours. 


Time  increment  in  decimal 

hours. 
Accumulated  decimal  hours.       SUM2 


Sy'mbo  1 
QSCALE  (I) 
SUM1 


Time  to  peak  in  decimal 
hours. 


Peak  discharge. 
Reciprocal  of  peak  discharge.     SSR 
Reciprocal  of  time  to  peak.  SSF 

Dimensionless  discharge  scale.  DEV 
Dimensionless  time  scale.  C 

Transformed  time  scale. 


SI  and  S2 
QC  (I) 


Descr  iption 
Transformed  discharge  scale. 
Uncorrected  sum  of  cross  products. 
Uncorrected  sum  of  squares. 
Parameters. 

Calculated  discharge  scale. 
Sums  of  squares  of  deviations  (rising). 
Sums  of  squares  of  deviations  (falling). 
Deviations,  observed  minus  calculated. 
Cut-off  constant. 
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OPTIONAL  SOURCE  PROGRAM  CHANGE 

The  Fortran  II  Source  Program  as  listed  in  Appendix  I  may  occasion- 
ally result  in  an  underflow  error  message  via  the  typewriter.   This  does 
not  affect  the  fitting  procedure  as  the  particular  hydrograph  ordinate, 
which  is  set  to  zero,  is  actually  near  zero,  i.e.,  EXP  (N)  where  N  <  -  227.9 
However,  one  may  program  around  this  by  making  the  following  source  program 
changes.   In  place  of  lines  102-103  inclusive  substitute  the  following: 

76  E  =  Sll-TSCALE(I) 
IF(E-f-200.)775,775,776 

775  QC(I)  =  0.0 
GO  TO  80 

776  QC(I)  =  EXPF(E) 
GO  TO  80 

77  E  =  S22*TSCALE(I) 
IF(E+200.)775,775,777 

777  'X(I)  -  EXPF(E) 


- 


APPENDIX  I 
FORTRAN  II  SOURCE  PROGRAM 


C      HYDROGRAPH  LAB  PROGRAM  NO.  1 

C      MATHEMATICAL  FIT  OF  A  DI MENSI ONLESS  HYDROGRAPH 

Q  ****************************************************************** 

Q  ****************************************************************** 

DIMENSION  TSCALE  (lOO).QSCALE  <  1  00  )  » DI SCH 1 < 100 ) .QC <  100 ) 
DIMENSION  TIME(  100) ♦ T I  ME  1(100)»TIME  2 (  100  )  »DI SCH ( 100 ) 
DIMENSION  51(10) .S2( 10) .  SSR<  10)  .SSF( 10) 

40  FORMAT  (  13. F3. 1.312. 14) 

41  FORMAT  (2F10.4) 
60  FORMAT  (315.15) 

62  FORMAT  ( 4F 10 . 4 » F 1 2 . 8 ) 

2  FORMAT  (2F2.0.I2.F7.4) 

3  FORMAT  (4F10.4) 

4  FORMAT  (2F10.4) 

Q  ******************************************************%•*********** 

25  READ  40.N.C.YEAR.M0.DAY.WTS 
q  ****************************************************************** 

C      CLOCK  TIME  TO  HOURS 
DO  103  1=1. N 
READ  2.T1.T2.T3»DISCH< 1  ) 
IF  (T3)998.6.5 

5  IF  ( Tl-12.  )8.999.8 
999  IF  (T2)6.8.6 

8  T1=T1+12. 
GO  TO  6 

998  IF  ( Tl-12.  )6. 997. 6 
997  IF  ( T2J996.6.996 
996  T1=T1-12. 

6  TIME! I )=T1+T2*. 016667 
IF  (  1-1  )101, 102.101 

Q  #»»*#■»*♦**♦#*****#**»*♦##**#***###*♦♦*♦*####******#*****#*»**♦#*** 

C  TIME    ACCUMULATION 

101  TIM=TIME( I )-TIME( 1-1) 
IF     JTIM17.7.995 

7  TIM=TIM+24. 
GO    TO    995 

102  TIMEK  I  )=0.0 
GO    TO    103 

995    TIME     MI)=TIME    KI-D+TIM 

103  CONTINUE 

Q  ****************************************************************** 

C      PEAK  AND  TIME  TO  PEAK 
TPEAK=TIME  1(1) 
PEAK=DISCH( 1 ) 
DO  10  1=2. N 
IF  (PEAK-DISCH(  I  )  )9.10.  10 

9  PEAK  =  DISCH( I  ) 
TPEAK=TIME  MI) 

10  CONTINUE 
C      **##«#*«*#*♦«*♦*******#»#***♦#*#*****#*****#**♦*******♦*#******#** 

C      OBSERVED  DI MENS  I ONLESS  HYDROGRAPH 

R=1./PEAK 

T=1./TPEAK 

DO  15  I=1.N 

DISCH  1  (  I  )=DISCH!  I ) *R 
15  TIME  2(1)  =  TIME  1  ( I ) *T 
C      **♦**##»♦#*#*#***#♦#******♦*»#*#****************************#***** 

C      TIME  AND  DISCHARGE  SCALE 
IF  ( TIME1 ( 1 ) 1400.400.401 

400  11=2 

GO  TO  601 

401  11=1 

11 
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601  IF  < DISCH1  IN)  > 500.500.50  1 

500  M  =  N-1 
GO  TO  402 

501  M  =  N 
402  DO  18  1=11. M 

A=LQGF( TIME  2(1)) 

B=SURTF(TIME  2( • ) ) 

B=l./B 

B=(B-1. )*2. 

TSCALE( I ) =A+B 
18  QSCALE( I )=LOGF(DISCH  1(1)) 
21  CONTINUE 

'*♦#************#*******#******************#*********************** 

PARAMETER  ESTIMATION 

KK1=1 

J1  =  0 
807  J1=J1+1 

SUM1=0. 

SUM2=0. 

DO  58  I  =  U»M 

IF  (QSCALEt I >+C)58.51 »51 
51  SUM1=TSCALE( I )*QSCALE( I )>SUM1 

SUM2  =  TSCALE(  I  >*TSCA|_El  I  )+SUM2 

IF( TSCALEf I )-. 000001 )  56»58»58 
56  Sll=SUMl/SUM2 

SUM1=0. 

SUM2=0. 
58  CONTINUE 

S22=SUM1/SUM2 

Sl( J1)=S11 

S2( J1)=S22 

CALCULATED  D I  MENS  I ONL ESS  HYDROGRAPH 
900  QC( 1)=0. 
J=l 

DO  80  1=11 «M 
IF( TSCALEl I )-. 000001 )  75.74.74 

75  J=J+1 
74  GO  TO  ( 76.77) .J 

76  QC(  I  )=EXPF(S11*TSCALEI 1  I  ) 
GO  TO  80 

77  QCt  I)*EXPF(S22*TSCALEU  >  ) 
80  CONTINUE 

IF(N-M)820.820.802 

802  QC(N)=0. 
820  GO  TO( 801. 16) »KK1 
801  SSR1=0. 

SSF1=0. 

DO  803  1=1. N 

DEV=DISCHl< I )-QC( I ) 

DEV=DEV«DEV 

SSR1=SSR1+DEV 

IF  (QC( I )-QC( 1+1 ) 1803.803.804 

803  CONTINUE 

804  DO  805  K=I »N 
DEV  =  DISCH1(K)-QC(IC) 
DEV=DEV*DEV 
SSF1=SSF1+DEV 

805  CONTINUE 
SSR( J1)=SSR1 
SSF( Jl )=SSF1 
C=C-.5 
IF(C)806.806.807 
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806  Jl=l  126 

JR=1  127 

XR=SSR(J1>  128 

D0  810J1=2»6  129 

IF<SSR( Jl >-XR)811»810»810  130 

811  XR=SSR(J1)  131 

JR=J1  132 

810  CONTINUE  133 

Jl=l  134 

JF=1  135 

XF=SSF(J1)  136 

DO  850  Jl=2,6  137 

IFISSFI Jl )-XF)851»850»850  138 

851  XF=SSF(J1)  139 

JF=J1  140 

850  CONTINUE  141 

PUNCH  72  142 

72  FORMAT  (1X.5H  YEAR.1X.3H  M0»1X»4H  DAY.2X»4H  WTS)  143 
PUNCH  60,YEAR»MO»DAY»WTS  144 
PUNCH  71  145 

71  FORMAT  (3X»5H  PEAK,3X»6H  TPEAK)  146 

PUNCH  41 »PEAK»TPEAK  147 

PUNCH  70  148 

70  FORMAT  (5X»5H  RISE»5X»5H  FALL)  149 

PUNCH  41.SKJR)  ,S2(JF»  150 

S11=S1(JR)  151 

S22=S2(JF)  152 

KK1  =  2  153 

GO  TO  900  154 

16  PUNCH  73  155 

73  FORMAT     <3X»6H    HOURS»5X,5H    T/TP»5X»2H    Q*6X.5H    Q/QP.5X.6H    OCALC )  156 
90    PUNCH    62.1 TIME1 I  I  )  »TIME2<  I ) »DISCH(  I  ) .DISCHK  I  )  ,OC«  I )  .  1  =  1, N)  157 

GO    TO    25  158 

END  159 
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APPENDIX  II 

INPUT  FORMATS 

HYDROGRAPE  LABORATORY  PROGRAM  NO.  1 

Parameter  Card 


Columns  Format  Symbol 

1  -     3  ^  XXXX  N 

4-6  XX.X  C 

7-8  XX  YEAR 

9-10  XX  MO 

11  -   12  XX  DAY 

13  -   16  XXXX  WTS 


Columns  Format  Symbol 

1  -     2  XX,  Tl 

3-4  XX.  T2 

5  -      6  +X  T3 

7-13  XXX.XXXX 


Description 
Number  of  lines  of  data  (max.  =  100) 
Cut-off  constant  for  QSCALE 
Last  two  digits  of  year 
Month  number  (1  -   12) 
Day  number  (l  -  31) 
Black  Book  watershed  number  (without  decimal) 


Data  Cards 

Hour  digit  of  time 

Minute  digit  of  time 

0,  military  time 
-1.  A.M.  and  noon 
+  1,  P.M.  and  midnight 

Discharge 


Description 


Note:    The  time  interval  between  discharge  tabulations  must  not  exceed  24  hours. 
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APPENDIX  III 


EXAMPLE  PROBLEMS 


001  is  parameter  card 


002  to  035  are  data  cards 


34  3  57  524  426 

— 

00  1 

251  1 

0056 

002 

252  1 

0144 

003 

256  1 

1300 

004 

258  1 

3420 

005 

301  1 

8300 

006 

304  1 

13000 

007 

307  1 

16800 

008 

310  1 

19400 

009 

312  1 

20800 

010 

315  1 

21800 

Oil 

317  1 

22000 

012 

319  1 

21700 

013 

321  1 

21300 

014 

324  1 

20200 

015 

327  1 

18500 

016 

330  1 

16500 

017 

336  1 

12500 

018 

340  1 

10000 

019 

345  1 

8240 

020 

355  1 

5400 

021 

401  1 

3530 

022 

405  1 

2990 

023 

410  1 

2680 

024 

422  1 

1840 

025 

439  1 

1150 

026 

455  1 

0756 

027 

515  1 

0509 

028 

650  1 

0144 

029 

800  1 

0072 

030 

1159  1 

0018 

031 

1200  1 

0018 

032 

100-1 

0015 

033 

200-1 

0010 

034 

300-1 

0005 

035 
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The  final  hydrograph  fit  QCALC  is  presented  for  the  rising  parameter  of  -3.7053  and  a  falling 
parameter  of  -8.1642.  The  actual  dimensionless  hydrograph  is  presented  in  the  4th  column. 


YEAR   MO 

DAY    WTS 

57     5 

24   426 

PEAK 

TPEAK 

2.2000 

•  4333 

RISE 

FALL 

-3.7053 

-8. 1642 

HOURS 

T/TP 

Q 

Q/QP       QCALC 

0.0000 

0.0000 

.0056 

.0025   0 

.00000000 

.0166 

.0384 

.0144 

.0065 

.00000001 

.0833 

.1923 

.1300 

.0590 

.03407449 

.1166 

.2692 

.3420 

.1554 

.  13408484 

•  1666 

.3845 

.8300 

.3772 

.3684854^ 

.2166 

.4999 

1.3000 

.5909 

.60572557 

.2666 

.6153 

1.6800 

.7636 

.78885133 

.3166 

.7307 

1.9400 

.8818 

.90837238 

.3499 

.8076 

2.0800 

.9454 

.95714221 

.3999 

.9230 

2. 1800 

.9909 

.99400277 

.4333 

.9999 

2.2000 

.9999   1 

.00000020 

.4666 

1.0769 

2. 1700 

.9863 

.98898750 

.4999 

1.1538 

2.1300 

.9681 

.9600011? 

.5499 

1.2692 

2.0200 

.9181 

,89444890 

.5999 

1.3846 

1.8500 

.8409 

.81472691 

.6499 

1.5000 

1.6500 

.7499 

.73051512 

.7499 

1.7308 

1.2500 

.5681    , 

.57021832 

.8166 

1.8846 

1.0000 

.4545 

.47729610 

.8999 

2.0769 

.8240 

.3745 

.37951267 

1.0666 

2.4616 

.5400 

.2454    . 

.23853391 

1.1666 

2.6923 

.3530 

•  1604     . 

.  18108410 

1.2333 

2.8461 

.2990 

.1359 

.15110291 

1.3166 

3.0385 

.2680 

.1218 

.12095815 

1*5166 

3.5000 

.1840 

.0836 

07226657 

1.7999 

4.1539 

.1150 

.0522 

03653790 

2.0666 

4.7693 

.0756 

.0343 

02018948 

2.3999 

5.5385 

.0  509 

.0231 

01020622 

3.9833 

9.1925 

.0144 

.0065 

00077044 

5. 1499 

11.8848 

.0072 

.0032 

00018107 

9.1333 

21.0774 

.0018 

.0008 

00000547 

9. 1499 

21.1158 

.0018 

.0008 

00000541 

10. 1499 

23.4236 

.0015 

•0006    . 

00000277 

11.1499 

25.7314 

.0010 

.0004 

00000150 

12.  1499 

28.0391 

.0005 

.0002 

00000085 
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